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Quantum gravity and spin systems* 
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Institut fiir Kernphysik, Technische Universitat Wien, A-1040 Vienna, Austria 

A new method for nonperturbative investigations of quantum gravity is presented in which the simplicial path 
integral is approximated by the partition function of a spin system. This facilitates analytical and numerical 
computations considerably. In two dimensions equivalence to an Ising model with ternary couplings is recovered. 
First simulations in four dimensions indicate strong similarities to the phase structure of original Regge theory. 



More than thirty years ago Regge developed 
a discrete description of General Relativity in 
which space-time is triangulated by a simplicial 
lattice, the Regge skeleton [1]. Thus, the lat- 
tice itself becomes a dynamical object, with the 
quadratic edge lengths q describing the evolution 
of space-time. Within this scheme the Einstein- 
Hilbert action for lattice gravity with cosmologi- 
cal constant A in c? dimensions is given by 

I(q) = X J2 V(s') - 2/3 ^ 6(.s''-^)V(.s''-'), (1) 

gd gd-2 

where the first sum is over all rf-simplices s** in 
the simplicial complex and V{s'^) is the rf-volume 
of the indicated simplex. The curvature of the 
lattice is concentrated on the (d — 2)-simplices 
leading to deficit angles 6(s'^~''^). j3 denotes the 
gravitational coupling. 

In the path integral formulation a quantization 
of the above action proceeds by evaluating the 
expression 

Z = j D[q]e-'^^) . (2) 

Unfortunately, a unique prescription for the mea- 
sure does not exist, however, an appropriate 
choice proposed in the literature is [2] 

/ = n/ 'fnui]), (3) 
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with m G J? defining a one-parameter family and 
T = \ for Euclidean configurations of squared 
edge lengths {</;} but T = Q otherwise. 

The central idea of this investigation is to trans- 
form the quantum gravity path integral to a par- 
tition function of a spin system [3,4]. The model 
is defined by allowing the squared edge lengths to 
take on only two values 

g; = H- ecT; , < e < emax , o-i e Z'j , (4) 

similar to the Regge-Ponzano model [5]. The real 
parameter e is restricted to fulfill the Euclidean 
triangle inequalities for the qi's so that all 2^^ 
configurations are allowed (Ni is the total number 
of edges), i. e. !F = 1 for all configurations {(/;}. 
Using (4) the measure (3) can be replaced by 

exp[-m ^ ln(l -|- ecr;)] = 
<7i=±i ; 

= ^ exp[-Vimo(e) - ^mi(e)(T;], (5) 
(7,=±i ; 

with niQ = — ime^-|-0(e'*) and mi = m(e+ ^e^) + 

In two dimensions, according to the Gauss- 
Bonnet theorem, the Einstein-Hilbert action as 
well as its discretized version is a topological in- 
variant. Therefore, it can be dropped if one con- 
siders manifolds with fixed topology. The only 
remnant of the gravitational action (1) is the cos- 
mological term. In the following the manifolds 
will be restricted to toroidal topology. The build- 
ing blocks of a 2-dimensional simplicial lattice are 
triangles and thus V(s'^) = V(s^) = At. To 
rewrite the action I = X Ai in terms of (T; 
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we consider a single triangle t with squared edge 
lengths qi,q2,qi. Its area can be expressed as 



At 



91 2(51+52-?;) 



(6) 



Expanding At the series consists only of terms up 
to (T^ since cr? = 1. This suggests the following 
ansatz 

At = Co(e) + Cl(e)((Tl + (72 + (T;) + C2(e)((Tl(T2 + 
+ (Tl(T; + (72(7;) + C3((f)(7l(72(7; . (7) 

There are only four possible values for the area 
of a triangle. Computing these areas and com- 
paring with (7) results in exact solutions for the 
coefficients c,-, e.g. 



C2 = -^[2^3- V(l-f)(3 + 5e) 
16 

-V(l + f)(3-5e)] . 



(8) 



Obviously one must have e < | = Cmax for the 
triangle areas to be real and positive. Inserting 
(5) and (7) into the partition function yields 



Z = J exp{— y~^(2Aci + mi)(7; — 

<7i=±i ; 

-A ^[C2((7l(72 + (71 (7; + (72(7;) + 

t 

-|-C3(7l(72(7;]}, (9) 

with J = exp(— A#2Co — A^i^o) and #2 the to- 
tal number of triangles. Thus, the path integral 
becomes the partition function of a system con- 
sisting of a spin (7; at each edge /, with an ex- 
ternal "magnetic field" and with 2- and 3-spin 
nearest neighbor interactions. Assigning the spin 
to the corresponding edge of the originally trian- 
gular lattice and drawing the interactions as lines, 
a Kagomis lattice is obtained. 

The term linear in a can be removed by a con- 
venient choice of the measure, mi = — 2Aci, and 
for e ^ the 3-spin couplings C3 can be neglected 
compared to C2. Hence, the remaining partition 



function is that of an Ising model on a Kagomis 
lattice 



Zr 



Jexp{(5 V'((7l-|-(72-|-(73-|-(74)(7;}, (10) 



<7l=±l 



where ai 2 and as 4 correspond to edges of neigh- 
boring triangles having (7; in common. To com- 
pute the critical coupling Qc = {—^C2X)c one 
uses the fact that Zi is connected to the par- 
tition function of a decorated honeycomb lattice 
via a triangle-star transformation and further to a 
conventional honeycomb lattice by a decoration- 
iteration transformation [6]. After all, one gets 

A. = -2^^, (11) 

C2 

which means that the antiferromagnetic regime, 
A > 0, is governed by frustration, whereas for 
ferromagnetic interaction, A < 0, the system un- 
dergoes a 2"''-order phase transition. Thus, close 
to the critical point our model with C3 = and 
nil = — 2Aci is equivalent to a (p^ field theory. 
On the other hand, 2d-gravity in the continuum 
is known to be equivalent to the Liouville theory. 
An analysis of our model (9) for reasonable mea- 
sure and couplings is necessary to shed light on 
the connection with continuum theory. 

A generalization of the 2-dimensional action to 
higher dimensions becomes rather complicated. 
In two dimensions only the three edges of a tri- 
angle are coupled in the action, whereas in three 
dimensions one is faced with terms of G*'' order at 
least, and with additional contributions from the 
Regge action. Recently, the Ising-link model in 
three dimensions has been analyzed via mean field 
techniques as well as with Monte Carlo methods 
[4]. However, critical exponents and the behavior 
of the curvature at the transition point /3c seem to 
differ from what was found for the Regge theory 
with continuously varying edge lengths. 

The situation is even more complicated in four 
dimensions where one has to deal with 10 edges 
in each simplex. Nevertheless, numerical simula- 
tions can be made very efficient using the heat 
bath algorithm. In the actual computations we 
replaced (4) by qi = &;(l-|-e(7;) allowing the lattice 
to fiuctuate around fiat space. A 2-dimensional 
simplicial lattice with equilateral links can be em- 
bedded in fiat space, but for the 4-dimensional 
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Figure 1. Spin expectation value as a function of 
l3 from simulations of spin quantum gravity on a 
toroidal S'* lattice. 
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Figure 2. Average Regge action for different cou- 
plings l3. Uniform measure m = 0, cosmological 
constant A = and e = 0.0875 are used. 



case hi has to take different values depending on 
the type of the edge /. In particular &; = 1, 2, 3, 4 
for edges, face diagonals, body diagonals, and 
the hyperbody diagonal of a hypercube, respec- 
tively. Preliminary studies show an interesting 
phase structure. 

The spin expectation value of a 4-dimensional 
toroidal lattice with S'* vertices is depicted in 
Fig. 1. The parameters were set to A = m = 0, 
e = 0.0875 and the range of /3 is extended to neg- 
ative values. The system undergoes phase tran- 
sitions at two different points, first in the negative 
region at /3 —4.67 and second at a large positive 
coupling /3 23.01. 

In Fig. 2 the expectation value of the corre- 
sponding Regge action as a function of /3 is dis- 
played. The resemblance of the curve for /3 > 
to that obtained from continuously varying edge 
lengths is very encouraging, cf. [2]. One has to 
determine the value of the critical points and the 
order of the phase transitions in the limit of in- 
finitely many vertices, Nq oo, in detail. Recent 
simulations indicate that both, the spin gravity 
model and Regge theory, belong to the same uni- 
versality class. If this is indeed the case one could 



perform calculations much easier than with con- 
ventional Regge calculus. 
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